Dynamics of the excitations of a quantum dot in a microcavity 
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We study the dynamics of a quantum dot embedded in a three-dimensional microcavity in the 
strong coupling regime in which the quantum dot exciton has an energy close to the frequency of 
a confined cavity mode. Under the continuous pumping of the system, confined electron and hole 
can recombine either by spontaneous emission through a leaky mode or by stimulated emission 
of a cavity mode that can escape from the cavity. The numerical integration of a master equation 
including all these effects gives the dynamics of the density matrix. By using the quantum regression 
theorem, we compute the first and second order coherence functions required to calculate the photon 
statistics and the spectrum of the emitted light. Our main result is the determination of a range 
of parameters in which a state of cavity modes with poissonian or sub-poissonian (non-classical) 
statistics can be built up within the microcavity. Depending on the relative values of pumping and 
rate of stimulated emission, either one or two peaks close to the excitation energy of the dot and/or 
to the natural frequency of the cavity are observed in the emission spectrum. The physics behind 
these results is discussed. 

PACS numbers: 78.67.Hc, 42.50.Ct 



I. INTRODUCTION 

Quantum electrodynamics of atoms within optical 
cavities is a well understood problem which has pro- 
duced many results of both fundamental and practical 
interest 0,13, 0H The current capability of using semi- 
conductors technology to grow quantum dots (QD) em- 
bedded in microcavities seems very promising to repro- 
duce and control that properties in a solid-state system 
that could be integrated in electronic or optical devices. 
The essential physics of such system is a strong cou- 
pling between the QD excitations and the photonic cavity 
modes as well as the possibility of interaction with the 
external world through the pumping of the system and 
the subsequent emission of light. The main goal is the 
control of the emitted light and its quantum properties 
being possible to built up high efficiency light emitting 
diodes, low threshold lasers and single photon sources. 

We concentrate here in a problem, no far from the ac- 
tual situation in several experiments [El Irj IsL Wi [Tol Hlf . 
in which a single QD is embedded in a pillar or disk 
microcavity confining photons in the three spatial direc- 
tions. The system is pumped either by electrically inject- 
ing electrons and holes which tunnel into the QD 6] , or by 
pumping excitons in the QD [E U EH EH ■ Light emis- 
sion form the system takes place both by cavity mode 
decay or by spontaneous emission of leaky modes. All 
these pumping and emission mechanisms introduce deco- 
herence affecting the quantum properties of the system. 

The lowest energy excitation of a neutral QD is an 
electron-hole pair, usually labelled as exciton. Due to 
the fermionic character of its components, this exciton 
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state can only be singly occupied having a degeneracy 
related to the possible values of the third component of 
the total angular momentum. This quantum number is 
usually referred to as the spin of the exciton. Among 
all these states, only two cases, those corresponding to 
±1, are of interest in our problem due to their possible 
coupling with photons. The next excitation corresponds 
to the case in which both +1 and —1 excitons are oc- 
cupied. This biexciton state has an energy different to 
twice that of a single exciton due to the Coulomb in- 
teraction of their components. This spectrum presents 
interesting alternatives ^2 that we intend to study in 
future work while here we restrict to simplest case, in 
which we consider only excitons with a given spin. This 
assumption is well justified in the case of experiments in 
which the system is pumped with polarized light, so that 
excitons with a given angular momentum are created in 
the system. We also assume that spin flip mechanisms 
are slow as compared with typical time scales in our sys- 
tem. Other single-particle states of the QD, above the 
exciton energy, can be neglected because they are not 
strongly coupled to the cavity mode in the case of res- 
onance or quasi-resonance between the exciton and the 
confined photon. Under these conditions, the problem 
reduces to that of a two-level artificial atom embedded 
in the optical cavity. We are going to show that, in spite 
of this simplification, the system presents a very rich va- 
riety of physical situations. 

The problem of a two-level system coupled to a single 
cavity mode, under different conditions and approxima- 
tions, has received a lot of attention, mainly in the field 
of quantum optics [JEHU ESQ 0f 

particular inter- 
est is the study of the possible sub-poissonian radiation 
when cavity losses and pumping dominate onto the spon- 
taneous emission of leaky modes. The aim of our work 
is an exhaustive analysis of the different regimes of pa- 
rameters to determine the role played by the different 
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physical mechanisms in the quantum properties of both 
the internal state of the system and the light emitted in 
the steady regime under continuous incoherent pumping. 

The paper is organized as follows: In section II we de- 
scribe the model Hamiltonian and master equation that 
allows to calculate the evolution of the populations and 
coherences of the energy levels. By means of the quan- 
tum regression theorem, we can use our master equation 
to calculate cavity photon statistics and the spectrum of 
the emitted light. In section III we present the results for 
the time evolution of the magnitudes of interest, in differ- 
ent ranges of parameters characterizing different physical 
situations. The spectrum of the light emitted by the sys- 
tem is presented and discussed in section IV. Section V 
is devoted to a summary of the work. 



II. THEORETICAL FRAMEWORK 

We consider a single QD inside a semiconductor mi- 
crocavity that is continuously and incoherently pumped. 
Our system is initially in its ground state and evolves 
until it reaches the steady-state. In order to describe the 
time evolution of the QD-cavity system we use a mas- 
ter equation that includes the strong exciton-photon cou- 
pling, the non-resonant pump, and the interaction with 
the environment that is responsible for dissipation. 



A. The Hamiltonian 

The physics of a QD strongly coupled with a single 
cavity mode is described by the following Hamiltonian: 

H = H S + H RS + H R (1) 

Hs, H R are the Hamiltonian for the QD-cavity system, 
and the environment (reservoirs) , respectively. The term 
Hrs describes the interaction between the QD-cavity 
system and the reservoirs. Hs is given by the usual 
Jaynes-Cummings Hamiltonian @> 3 that describes the 
interaction of a two-level system with a single mode of 
the electromagnetic field: 

Hs = H o + Hx-c 
H Q = co x \X)(X\ + (Lu x -A)a^a 
Hx-c = g{cra)+a^) (2) 

We have introduced ladder operators <r =| X)(G \, 
a =| G)(X | connecting the ground \G) and excited (ex- 
citon) \X) states of the QD with energies zero and u>x 
respectively (we take h = 1). creates a cavity photon 
with energy (ojx — A). The coupling term Hx-c in Eq. 
@ describes the exciton-photon coupling in the Rotating 
Wave Approximation 0, 0, . 

Hrs contains the coupling to external reservoirs in- 
cluding the following three process: 



(i) The continuous and incoherent pumping of the QD 
by annihilating (c R >) an electron- hole pair from an 
external reservoir R' (representing cither electrical 
injection or the capture of excitons optically cre- 
ated at frequencies larger than the typical ones of 
our system) and creation (d RII ) a phonon emitted 
to a reservoir R" in order to take care of energy 
conservation 

^ ^w,R"[d R ,,c Rl ^ + ac R ,d,R»], (3) 

R',R" 

(ii) The direct coupling of the QD exciton to the leaky 
modes, that is, to the photonic modes, with energy 
different than the cavity mode, that have a resid- 
ual density of states inside the microcavity. This 
process is responsible for the dissipation of the ex- 
citonic degrees of freedom by the spontaneous emis- 
sion to an external reservoir R of photons created 

by^ 

J2*lR[°b R + b R ai], (4) 

R 

(iii) The escape of the cavity mode out of the microcav- 
ity due to the incomplete reflectance of the mirrors. 
The cavity mode is thus coupled to the continuum 
of photonic modes out of the microcavity. This 
process produces the direct dissipation of the cav- 
ity mode 

H CE = Y, X R( ab R + b R a1 )- ( 5 ) 

R 

The last term in Eq. H R , describes the exter- 

nal reservoirs of harmonic oscillators (photons, phonons, 
electron- hole pairs,...) not being necessary to detail them 
explicitly The three terms H R s have coupling constants 
I^R',R", Ai? and Xi R which depend on the particular mode 
(R, R' or R") of each external reservoir. The operators 
&j b R , c R i and d R n have bosonic commutation rules as 
it corresponds to harmonic oscillators. Since we are in- 
terested in the strong coupling regime, the first three 
terms, Hq and Hx-c ar e treated exactly. That means 
that one could work in the framework usually known as 
the "dressed atom picture". However, in order to clar- 
ify the different pumping and losses mechanisms, it is 
preferable to work in basis {| Gn); | Xn}} in terms of 
the number of cavity photons n and bare ground G and 
exciton X states of the QD. 

B. Master equation 

We have made the whole algebra for a general case in 
which reservoirs are at finite temperature. However, the 
main physics already occurs for zero temperature, which 
is the case we present hereafter in this paper. 
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We define the reduced density matrix, p, for the 
exciton-photon system by tracing out the reservoir de- 
grees of freedom in the total density matrix px- 



Tr R (p T ). 



(0) 



In the interaction picture with respect to Hx-c + Hrs 
in Eqs. and J3J, p satisfies the master equation 0, Ell : 



d_ 

dt' 



[p,H s } + 
— (2apa) — a)ap — pa) a) - 
(2apa^ — a' up — pa' a) 
(2a* pa — aa^ p — paa') 



(7) 



The master equation is obtained under the usual Born- 
Markov approximation for the interaction Has between 
the QD-cavity system and the reservoirs, but the strong 
exciton-cavity photon coupling, Hx-Ci is described ex- 
actly, k is the decay of the cavity photon by escaping 
through the microcavity mirrors, 7 is the decay of the QD 
exciton by the spontaneous emission into leaky modes, 
and P is the rate of continuous incoherent pumping of 
the QD exciton. By means of Eq. one can get a 
set of differential equations that describe the evolution 
of the populations and coherences of the cavity-QD sys- 
tem. In the basis {\Gn)\ \Xn)} of product states between 
QD states and Fock states of the cavity mode, the matrix 
elements of the reduced density matrix are: 



Pin,jm 



(in\p\ji 



(8) 



with i,j being either G or X. The diagonal matrix el- 
ements PGn,Gn, PXn.Xn are the populations of the QD- 



photon levels, while the non-diagonal terms PGn,Xn-i, 
PXn-i,Gn, describe the coherences between these levels. 
By taking the matrix elements in Eq. J7J we get, for 
T = 0, the following set of linear differential equations: 
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FIG. 1: Ladder of levels (continuous lines) for a two-state 
QD coupled to a single optical mode of a microcavity. States 
labeled as Gn with n = 0, 1,2, ... correspond to having the 
QD in its GS coexisting with n cavity photons. The same for 
states Xn with the QD in its excited state X. The double, 
dashed and dash-dotted lines depict the coupling, pumping 
and emission processes with rates g, P, 7 and k explained in 
the text. 



dtPGn,Gn = igV™ [pGn,Xn-l ~ PXn-l,Gn) 

+"IPXruXn - K [npGn,Gn ~ {n + l)pGn+t,Gn+l] - PpGn,Gn (9) 

dtPXn,Xn = igVn+ 1 (pXn,Gn+l ~ PGn+l,Xn) 

-JPXruXn - « [npXn,Xn ~ (« + 1) PXn+l,Xn+l] + PpGn,Gn (10) 

dtPGn,Xn-l = i [gVn (PG - PXn-l,Xn-l) + ApGn,Xn-l] 

- [(7 + K(2n - 1) + P) /2] PGn,Xn-l + K\/n{n+l)pGn+l,Xn, (H) 

plus the equation hermitian conjugate of In the absence of dissipation, the matrix elements PGn.Gn, PXn-i.Xn-i, 
PGn.Xn-i, and pxn-i.Gm for a given photon number n, satisfy a closed set of four differential equations. However, 
the pumping and emission with rates P, k and 7, couple the terms with different photon occupation number n, so 
that an infinite set of equations has to be solved, as depicted schematically in Fig. 1. For the numerical integration, 
the set of equations can be truncated at a given value that, in all the cases considered below, it is enough to take 
equal to 100. 

As initial conditions we take pco,Go = 1 and all the other elements of the density matrix equal to zero. In other 
words, we start with the system in its ground state and the pumping, and subsequently the losses, produces the 
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dynamical evolution of the whole system. The steady-state of the system can be studied by integrating the set of 
equations ©, I|1U|) and Ijllfl for long times. 



C. First and second order coherence functions 



From the master equation, one can compute the dynamics of the expectation values of any operator. Moreover, 
two-time correlation functions are also of physical interest. In particular we want to compute the first and second 
order coherence functions 



y/(E(~) (f, t)E(+) (r, t)) (£F) (r, t + t)EW (r, t + r)) 

(2) (gM (r , t)E(-) (r, t + r)£(+) (r , t + t)£<+> (r , t)) 

5 1 ' ' J (£(-) (r, *)£(+) (r,i))(£(-)(r,i + T)£(+)(f,t + r)) [ ' 



where E^ + ' and £A _ ) are the positive and negative fre- 
quency parts of the amplitude of the electromagnetic 
field. In the steady-state limit we are interested in, both 
and do not depend on the absolute time t. 
The experimental situation is such that photonic 
modes escaping from the cavity are emitted in well de- 
fined direction, for instance along the axis of the micropil- 
lars in the direction in which the mirrors defining the 
cavity have some 

transparencyBHIlllBSIijI. On 

the contrary, emission of leaky modes takes place in any 
direction, particularly through lateral surfaces of the mi- 
cropillar. This means that simply by changing the spatial 
distribution of detectors, one can measure the emission 
from cavity modes or the leaky modes. 

By Fourier transforming the first order correlation 
function 

G«(t,r) = {E^(f,t)EW(f,t + r)), (14) 

one can obtain the the power spectrum of the emitted 
light. 

1 f°° 

S{f,v) = -fft dre il,T G w (t,T). (15) 

I" JO 

The first order correlation function of the external field 
can be obtained from the time dependence of the opera- 
tors describing intrinsic properties of the system: 

G^(t,r) = (aHt + r)a(t)) 

G«(t,r) = (^(t + r)a(t)). (16) 

G(p (t, t) is the correlation function for the cavity mode, 
responsible for the stimulated emission part of the spec- 
trum (i.e. the light coming from the confined photon) 

by means of G^\t,T) cx kGq (t,r). In the case of 
pillar microcavities, it would correspond to the light 
emitted in the vertical direction. On the other hand 
G)P {t, t) describes the spontaneous part of the spectrum 



(i.e. light directly coupled to the QD exciton) by means 

of G {1 \t,r) cx 7 G^ ) (^ T )- This gives the light emitted 
through the leaky modes that can be measured in the 
lateral direction of a micropillar. 

In order to calculate these two-time correlation func- 
tions, we make use of the quantum regression theorem 
First of all, we define the following operators: 

4» = \Gn + l){Gn\e^-^ 

a Xn (r) = \Xn + l)(Xn\e^-^ 

4(r) = \Xn)(Gn\e^ 

fc,(r) = \Gn + l)(Xn-l\e^- 2A ^. (17) 

In the interaction picture, time evolution with respect 
to Hx-c + HrS appears explicitly. The two-time cor- 
relation function of the cavity mode can be expressed in 
terms of these operators: 

G { c\t,T) =^V^TTx (18) 

n 

(<ok(* + *>(*)) + ( a xJt + r)o(t)>) ■ 
The exciton correlation function can be written as: 

G^(t : r)=J2(^(t + r)cr(t)) (19) 

n 

For the calculation of the spectrum of the emitted light 
it is, thus, necessary to evaluate the functions (a>Q n (t + 
r)a(t)), (a Xn (t + T)a(t)), and (<r+ (t + T)a(t)). The quan- 
tum regression theorem states that given a set of oper- 
ators Oj, whose averages satisfy a closed set of linear 
differential equations: 

±{ 0j (t + T)) =J2 L iMOk(t + T)), (20) 
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then the two-time averages of Oj with any other operator 
O, also satisfy the same differential equation: 

±( 0j (t + T)0(t)) =Y,h,k{Ok{t + r)0{t)) (21) 
fc 

In order to get the time evolution of two-time averages, 



we start with the dynamics of the operators l|17(l . Their 
averages satisfy the set of linear differential equations 
that allows us to find the evolution of the corresponding 
two-time averages in <|19l I19f) : 



8r<4»W> = (dr P Gn,Gn + l) e^~^ + i(« - A)<4„(t)> 

Sr<4W) = {drPGn,Xn)e^ T +lUj(al{T)) 

dA^Xn-Ar)) = (8rPxn-i,*n) e^ A ^ + i(w - A)(4„_ 1 (r)) 

5 r (ft,(T)) = (d T p Xn - hGn + 1 )e^- 2A ^ +i(w-2A)(ft,(T)). 



(22) 



Although the operator <j n is not needed in the evaluation of the required two-time averages (|19|) and 119( 1 . we have to 
add the correlation functions that include such operator <j„ in order to get a closed set of equations. Using the master 
equation and eliminating the elements of the density matrix, we arrive at the desired time-evolution for the averages 

dr{a Gn {r)) = [i("x - A) - |(2n + l) - P 

+ W(" + + 2) + (4(r))*5>/Jm 

+( fl L-i( T ))7- (s n (T))igVn 



d T (al(r)) = (a Gn (T)^gV^+T+ (4(r)> 



(7 + P) 

ZCiJx KU 



+(<4+i( r )) K ( n + !) - (4cn-i( T )) i 5V™ 

+<«L-r M) - A) - 7 - ~(2n - 1) 
+ < a x„( T )) K vM™ + 1 ) + (cW)iffVn+l 

9 T (?„(r)) = -( a ( G „)( T )) J sVn + (a{ Xn _ 1) (T))ij\ArrT 



(?n(r)) 



n w x - A) nil 



-{^ n +i(r))Ky/ n(n + 2). 



(23) 



From Eq. (|23|l . and the quantum regression theorem, it is straightforward to obtain two separate closed sets of 
differential equations for two-times functions: one for the set of functions 



{(aln(r + *)<*(*)>> («k-i(^ + t)a(t)), (4(r + t)a(t)), ( fe (r + t)a(f)»> 
needed for the calculation of G^\t,r), and other for the set of functions 

{(aljr + t)a(t)), {a} Xn _ x {r + t)a(t)), (4 ( T + t)^)), ( ?n ( T + t)<r(t)>}, 



(24) 



(25) 



needed for the calculation of G x (t, r) 



An important point is the initial conditions to solve these systems of equations. Such conditions are obtained by 
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solving, up to the stationary limit, the master equations 
© , (imp and for the density matrix. From this infor- 
mation, one knows the initial conditions. For the func- 
tions in (J23J: 

( a Gn(*) a (0> = Vn + lpGn+l,Gn+l 
( a Jf«-l( i ) a ( i )) = VnpXn,Xn 

(o-t(t)o(t)) = Vn + lpGn+i,x n e~ iM 



iAt 



(26) 



For the functions in 1251) : 



( a Gn(*) <7 (*)) = PXn,Gn+lt 

(ok-i(*M*)> = 

(&,<7(t)) = 



-iAt 



(27) 



The second order coherence function is, a priori, more 
complicated to calculate. Averages of products of four 
operators at two different times have to be performed. 
This task becomes much simpler for the case of zero time 
delay. g( 2 >(t,T — 0) is a one-time operator, which sim- 
plifies the calculation. In spite of loosing information, 
g^(t,r = 0) is a very interesting magnitude because it 
can be used as indicator of the possible coherence of the 
state of the systempj ESQ. In addition, if we concen- 
trate in the properties of cavity modes (i.e., neglecting 
the emission of leaky modes), the second order coherence 
function acquires a very simple form in the stationary 
limit: 



3 (2) (i,r = 0) 



(a'a'aa) 

(flt a )2 

J2n n ( n ~ l )[PXn,Xn + PGn,Gn 
En n [PXn,Xn + PGn,Gn]] 2 



(28) 



<?( 2 ) takes different values depending on the statistics of 
the photon state: = 2 for chaotic states, g^ = 1 for 
states having a Poisson distribution in n, and g( 2 ' < 1 for 
non classical systems having a sub-Poisson distribution 

aaa. 



III. RESULTS 
A. Election of the range of parameters of interest 

Since our model has several parameters, we need some 
simple picture in order to get insight on the interesting 
regime of parameters. Although our system is similar to 
a spin coupled to an harmonic oscillator, we can under- 
stand the effect of dissipation by considering two-coupled 
(by g) oscillators. One of them has an eigenfrequency 
ujx — A being damped with a rate k while the other, 
with an eigenfrequency cox, is damped with a rate 7 
and is pumped incoherently with a rate P. The inco- 
herent pumping implies a dephasing mechanism which, 



together with the damping, gives an average dephasing 
(7 + P)/2. When the damping of the first oscillator k 
is much higher than the dephasing (7 + P)/2 of the sec- 
ond oscillator, one could expect that the frequency of 
this second oscillator lox will survive longer than ujx — A 
and it will dominate the spectrum. In the opposite limit, 
k <C (7 + P)/2, the frequency lox ~ A will dominate 
the spectrum. This behavior is the one shown by our 
spin-oscillator system as can be checked by means of the 
Haken's model on adiabatic elimination 0, 01 0] in the 
equations of motion of the field operators. From this 
simple picture we can motivate the existence of three dif- 
ferent regimes n ^> (7 + P)/2, n -C (7 + P)/2 and some 
intermediate regime. For this purpose, we fix g = 1 as 
the energy scale, take a fixed small value 7 = 0.1 reflect- 
ing the fact that leaky modes emission is the less efficient 
mechanism in practical situations, and define the three 
regimes by changing k and P. 

Apart from those parameters, our model requires to set 
the frequency of the two oscillators. In practical cases, 
the coupling between the oscillators is in the order of 
meV while the excitation energy lux is in the order of 
eV. Therefore, since our energy scale is g — l(meV), all 
the calculations we present here have been performed by 
taking lox = 1000. Finally the detuning A is a very 
important parameter. Since typical detunings are in the 
order of a few meV, we present here results for two differ- 
ent cases: perfect resonance, A = and quasi-resonance 
A = 5. 



B. Dynamics of the density matrix: occupations 
and coherences 



In this subsection we show the time evolution of the 
occupations and coherences (diagonal and off-diagonal 
respectively) Pinjm described in Eqs. (P|). l(TU|> and (fTTfl . 
From these equations, it is clear that if losses and pump- 
ing are not included, the systems shows the usual Rabi 
oscillations in each subspace with a given number of pho- 
tons. The amplitude of these oscillations increases with 
the coupling, <?, and decreases with the detuning A. 

When all the losses and the pumping are included, the 
situation changes drastically. The initial occupation of 
I GO) evolves with time up to a state of equilibrium in 
which occupation is redistributed among a large number 
of levels with a finite number of photons. In this final 
steady situation, the sum of all the losses equilibrates 
the pumping. 

Figs. 2 and 3 show the time evolution of occupa- 
tions pin,in for both the exciton (i = X) and the ground 
(i = G) states. The results in Figs. 2 and 3 are typi- 
cal of the two extreme regimes described above. When 
pumping rate P dominates (Fig. 2), states with a large 
number of photons n become occupied and two bundles 
of branches can be distinguished: upper branches corre- 
sponding to I Xn) states and lower branches correspond- 
ing to I Gn) states. In any case, the differences between 
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FIG. 2: Density matrix diagonal elements pi n ,in for both the 
exciton (i = X) and the ground (i = G) states, from n — 
up to n = 35 cavity photons, g = 1, A = 5, 7 = 0.1, K = 0.1 
and P = 15. The inset shows a zoom of the graphic to clarify 
the formation of two bundle of branches associated with the 
exciton and the ground states of the QD respectively. 



the occupations of states in the two bundles, is not very 
large. When rates for losses increase and become com- 
parable to the pumping rate, just a few photons can be 
stored in the cavity. In equilibrium, the highest occu- 
pation is that of the | GO) state and only states with 
n < 3 have appreciable, although small, occupations. It 
is quite evident that an interesting regime is missing in 
these results: the one in which losses dominate onto the 
pumping, for instance the case k = 5 and P = 1. Wc 
will discuss below other characteristics of this case but 
do not include a figure here because the corresponding 
time dependence is featureless showing a fast decay to 
zero of the occupations for all the states with n > 0. 

Time evolution of coherences pxo,Gi is shown in Fig. 
4 for the two limiting cases P ^> k and P <C k for A = 5. 
Due to this finite detuning, Rabi-like oscillations appear 
in the coherences. When k dominates, the coherences in- 
crease with time up to finite, but small, value. When 
pumping dominates, the mean number of photons in- 
side the cavity is much higher than 1 as discussed below. 
Therefore, a coherence like the one in Fig. 4 between 
states with n = and n = 1 decreases with time. In this 
case, P ^> K, we have computed coherences pxn,Gn+i, 
for values of n around the mean number of photons in- 
side the cavity (see below), obtaining that they also go 
to a small constant (apart from the Rabi oscillations) for 
large times. 

All the results included in Figs. 2, 3 and 4 correspond 
to a finite detuning A = 5. In the case of perfect reso- 
nance A = 0, the main effect is the lack of oscillations 




t (psec) 



FIG. 3: Density matrix diagonal elements pi n ,in for both the 
exciton (i = X) and the ground (i = G) states, from n = 
up to n = 3 cavity photons, g = 1, A = 5, 7 = 0.1, k = 0.5 
and P = 1. The inset shows a zoom of the vertical axis of the 
graphic. 
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FIG. 4: Density matrix off diagonal element pxo,Gi with de- 
tunning A = 5 for g = 1, 7 = 0.1, and two different values of 
k and P. 

as shown, as an example, in Fig. 5 for the coherences 
Pxo,gi- 

The fact, shown above, that P competes with k, is not 
completely general. This fact can be seen in equations 
ijffil. (|10|l and (flip for the time evolution of the density 
matrix elements. For the off-diagonal terms, P has the 
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FIG. 5: Imaginary part of the density matrix off diagonal 
element pxo,Gi at perfect resonance (A = 0) for g = 1, 7 = 
0.1, and different values of k and P. All the real parts of the 
same matrix elements are zero. The inset shows a zoom of 
the graphic for small values of t. 



same sign than the losses k and 7, while for the diagonal 
terms, these signs are different. The pumping produces 
decoherence in the off-diagonal terms as the emission of 
photons does because in our model the pumping is in- 
coherent. At the same time, P favors higher occupa- 
tions (increase of the diagonal elements). Therefore it is 
clear that, for a high pumping, decays of the off-diagonal 
terms (loss of coherence) prevail on the occupation until 
the coherence is completely quenched and later inverted, 
preventing the evolution towards the occupation of states 
with a high number of photons. This is an indication of 
the important role that off-diagonal elements play in the 
time evolution. It also explains that, if the emission rate 
k is above a critical value, it dominates on any pumping 
effect. 



C. Number of cavity photons 

A very important result to be drawn from the time 
evolution discussed above is the fact that the system is 
able to store a significant number of cavity photons. Let 
us analyze the mean number of cavity photons as a first 
step to study the quantum properties of these photons by 
looking to its distribution. The mean number of photons 
in the cavity is 

N P h = X! n (PXn,Xn + PGn.Gn)- (29) 
n 
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FIG. 6: N p h, in grey scale, as a function of n and P for perfect 
resonance A = with g = 1 and 7 = 0.1. 
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FIG. 7: N p h, in grey scale, as a function of k and P for a 
detuning A = 5 with g = 1 and 7 = 0.1. 

N p h increases with time, in a scale of hundreds of pi- 
coseconds, up to the stationary value that is the magni- 
tude we are interested in. Figs. 6 and 7 show a contour 
plot of the mean number of cavity photons for A = 
and A = 5 respectively. As in previous figures, we fix 
all the parameters except P and k, which vary along the 
two axis. In these figures a region with a high N p ^ can 
be observed. It corresponds to low emission and high 
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FIG. 8: Second order coherence function g^ 2 \r = 0), in grey 
scale, as function of k and P for perfect resonance A = 0. g = 
1 and 7 = 0.1. The dashed area shows the region supporting 
a Poisson distribution of cavity photons in which <T 2 '(t = 0) 
close to 1. 



FIG. 9: Second order coherence function g^(r = 0), in grey 
scale, as function of k and P for a detuning A = 5. g = 1 
and 7 = 0.1. The dashed area shows the region supporting 
a Poisson distribution of cavity photons in which <T 2 '(t = 0) 
close to 1. 



pumping. The region is larger and the number of cav- 
ity photons higher for the resonant case, A = 0, than 
for the case with detuning A = 5. A detuning makes 
more difficult the coupling between the states | Xn) and 
| Gn +1). As a consequence the mechanism of storing 
cavity photons become less efficient. 

A conclusion similar to the one drawn in the previous 
section is in order: when pumping increases, the mean 
number of cavity modes also increases up to a certain 
maximum value. Further increase of P implies that the 
dephasing mechanism related with this incoherent pump- 
ing provokes a decrease of N p h- 



D. Second order coherence function 

The possibility of storing a large number of cavity pho- 
tons for a wide range of parameters opens an interesting 
alternative beyond the simple value of N p h' how is the 
distribution of cavity photons? If any interesting distri- 
bution is involved, its character could be transfered to 
the light emitted by the system. 

In order to quantify if such distributions are closer to 
gaussian distributions (thermal or chaotic states), Pois- 
son distributions or even non-classical sub-poissonian dis- 
tributions, one can compute, using Eq. (|29[) . the second 
order coherence function g^ 2 \r = 0) discussed above. 

Figs. 8 and 9 show ,g (2) (0), for A = and A = 5 
respectively. In these two figures we have dashed the 
region in which g^ 2 \0) is close to 1 to show the region 



supporting a poissonian distribution. This is the border 
between classically accessible region and the one having 
gW (0) < 1 with non-classical states with sub-poissonian 
distribution. In the case of perfect resonance (Fig. 8), for 
low k, the poissonian distribution region is rather wide in 
terms of the pumping P. When n increases, this region 
becomes much narrower in P. In the case of detuning 
A = 5, the poissonian region becomes much smaller and 
even disappearing for intermediate values of k. 

The interesting region of g^ (0) < 1 appears for very 
low pumping in these two figures. Even though this oc- 
curs in a small region of parameters, it is very interest- 
ing because it means that the system can form squeezed 
photon states with sub-poissonian distribution. In this 
region the emitted light shows the antibunching that is 
characteristic of non-classical light emitters. 

The rest of the diagram is supporting states with 
(/ 2 )(0) increasing from 1 to 2, i.e. states where second 
order coherence is reduced and approaching a gaussian 
distribution (g^ 2 \0) = 2). A finite detuning A enlarges 
this gaussian region as observed when comparing the two 
figures. 

Combining and summarizing the results of Figs. 6, 7, 8 
and 9, it must be stressed that, for high quality cavities, 
i.e. low 7 and k, and intermediate value of the incoher- 
ent pumping P is able to produce a rather large num- 
ber of cavity photons with Poisson distribution. When 
detuning increases this effect remains, although with a 
reduction of the region of parameters supporting it. 
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FIG. 10: Envelope of the first order correlation functions 
(t) (upper panel) and G^\t) (two lower panels) at the 
stationary limit (t — > oo). g = 1, A = 5, 7 = 0.1, k = 0.1 and 
P = 15. 



FIG. 11: Envelope of the first order correlation functions 
G^(t) (upper panel) and G^\t) (two lower panels) at the 
stationary limit (t — > 00). g = 1, A = 5, 7 = 0.1, k = 5 and 
P= 1. 



IV. SPECTRUM OF THE EMITTED LIGHT 

In this section we will present our results for the spec- 
trum of the emitted light in the stationary limit (t — ► 00). 
This is the property that can be most easily measured 
in experiments. We will do it only in the case of finite 
detuning A = 5 in order to observe the emission at two 
different frequencies. Following the steps discussed in the 
section devoted to our model, we have computed G^\t) 
and Gq^ (t) in the stationary limit (i — > 00) for the three 
regimes already discussed P> k, P <C k and P com- 
parable to k. The envelopes (with fast oscillations in- 
side) of these functions are shown in Figs. 10, 11 and 
12 in which we show different contributions to the to- 
tal emission. In these three figures, the upper part shows 
the spontaneous emission of leaky modes (G^ (t)), while 
the two lower parts show two different contributions to 
the stimulated emission through cavity modes (G^'(r)). 
Fig. 10 correspond to a case with a large N p h inside the 
cavity. Therefore, when summing up in n to get G^ (r) 

and G { c\t), one obtains a time scale significantly larger 
than in the case of Figs. 11 and 12, which correspond to 
a small N p h inside the cavity. 

The Fourier transform of Figs. 10, 11 and 12 gives the 
spectrum of the light emission in the three cases discussed 
above. In order to numerically perform such Fourier- 
transforms, we have used digital data filters (Parzen and 
Haning) to reduce numerical noise. The spectra corre- 
sponding to lox = 1000 are shown in Figs. 13, 14 and 
15 . In these three figures we have made use of the fact, 



discussed in the introduction, that the two types of emis- 
sions can be separated: the emission coming from the 
cavity modes takes place along the axis of the pillars, 
while the emission of leaky modes takes place in any spa- 
tial direction (mainly in the directions perpendicular to 
the axis of the pillar). 




FIG. 12: Envelope of the first order correlation functions 
Gjc( t ) ( u PP er panel) and G^?(t) (two lower panels) at the 
stationary limit (t — ► 00). g = 1, A = 5, 7 = 0.1, k = 5 and 
P= 15. 
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FIG. 13: Spectrum of emission for g = 1, u>x = 1000, A = 5, FIG. 14: Spectrum of emission for g = 1, u>x = 1000, A = 5, 
7 = 0.1, k = 0.1 and P = 15. 7 = 0.1, k = 5 and P = 1. 



The case of Fig. 13 corresponds to a high pumping 
dominating on the effect of k. As discussed above, one 
expects a strong damping of the oscillator with frequency 
lux- Therefore, only one peak at lox — A is observed. 
Moreover, the peak shown in this spectrum is signifi- 
cantly narrower than the features observed in any other 
case we have analyzed as, for instance, the ones in the 
following figures. 

Fig. 14 shows the result with a high value of the rate 
emission n — 5 and a smaller, although non negligible, 
pumping P = 1. This case corresponds to N p h = 0.027 
and </ 2 -*(0) = 0.396, i.e. to a sub-poissonian distribution 
of a small number of cavity photons. Now, the strong 
damping of the mode with frequency u>x — A reduces its 
intensity which becomes smaller than a second peak at 

LU X - 

In the intermediate case shown in Fig. 15, the dephas- 
ing of the two modes is similar and this produces a very 
wide peak at the spectrum as observed in the figure. 

Let us finally analyze more carefully the case of high 
quality cavity (already studied in Fig. 13). For this pur- 
pose, we maintain fixed a low value k = 0.1 and let the 
pumping vary in a wide range of values. Figs. 16 and 
17 show the spectra of the stimulated (cavity) and spon- 
taneous (leaky) emissions respectively. First of all one 
must observe that scale in Fig. 16 is more than 200 times 
larger than that of Fig. 17 as it could be expected from 
the results and discussion of Fig. 13. The main result to 
be drawn from Figs. 16 and 17 is that the most intense 
emission corresponds to the range 10 < P < 30 simply 
because this is the range of higher number of photons 
inside the cavity as shown in Fig. 7. 
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FIG. 15: Spectrum of emission for g 
7 = 0.1, k = 5 and P = 15. 



1, lox = 1000, A = 5, 



V. SUMMARY 

We have described the dynamics of a QD embedded 
in a semiconductor microcavity by means of a two-level 
system strongly coupled to a single cavity mode. The 
system is continuously excited by incoherent pumping of 
excitons. Two different sources of photon emission are 
considered: spontaneous emission through a leaky mode 
and stimulated emission of a cavity mode escaping from 
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FIG. 16: Spectrum, in grey scale, of stimulated emission (cav- 
ity modes escaping from the cavity) as a function of frequency 
and P. g = 1, u x = 1000, A = 5, 7 = 0.1 and k = 0.1. 



served. When the rate emission k is higher than a non 
negligible pumping rate P, the strong damping of the 
mode with frequency lox — A allows to observe the two 
modes with a higher intensity for the mode at iox- In 
a intermediate case, the dephasing of the two modes is 
similar producing a very wide peak at the spectrum. 




frequency (meV) 



the cavity. The time evolution of first and second or- 
der coherence functions is calculated. When pumping 
dominates over emission rates, a large number of cavity 
photons can be stored in the cavity. Further increase 
of the pumping introduces dephasing and a decrease of 
the number of cavity photons. These different regimes 
are also characterized by poissonian or gaussian photon 
distributions inside the cavity. Sub-poissonian distribu- 
tions can be obtained for a range of parameters, in which 
the pumping rate is very small, and the quantum nature 
of the QD-cavity system manifests itself in the emitted 
light. Finally, we have studied the emission spectrum of 
our system. In the case of high pumping dominating on 
the rate emission k, one gets a strong damping of the os- 
cillator corresponding to the exciton level with frequency 
u>Xi an d only one very narrow peak at lux — A is ob- 



FIG. 17: Spectrum, in grey scale, of spontaneous emission 
(emission of leaky modes) as a function of frequency and P. 
g = l,ux = 1000, A = 5, 7 = 0.1 and k = 0.1. 
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